Oscillatory nonlinear conductance of an interacting quantum wire with an impurity 
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The nonlinear conductance of a one- dimensional quantum wire adiabatically coupled to Fermi 
liquid electron reservoirs is determined in presence of an impurity. We show that electron-electron 
interaction in connection with the finite length of the wire leads to characteristic oscillations in the 
current as a function of the applied voltage. 
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It is well known that Fermi Liquid (FL) theory is 
not appropriate to describe transport properties of one- 
dimensional (ID) electron systems. Correlation effects 
completely destroy Landau's quasiparticle picture, and 
the elementary excitations of a ID quantumwire (QW) 
can be modelled in terms of bosonic plasmon and spin 
density modes. An equivalent representation employs 
fractionally charged quasiparticles |ll The charge 
ge depends on the parameter g < 1 characterizing the 
electron-electron interaction strength, where 5 = 1 cor- 
responds to the non-interacting limit and g means 
strong correlations. An impurity potential in a ID QW 
is strongly renormalized by the interaction 3] and ef- 
fectively acts as a back-scatterer of these fractionally 
charged quasiparticles Ij. 

These well known properties of ID fermions are a natu- 
ral consequence of the Luttinger Liquid (LL) model ^|^, 
which generically describes the low energy physics of 
QWs. Despite these drastic departures from FL theory, 
the strong correlation effects are often disguised in trans- 
port measurements. For instance, the conductance of a 
clean QW is not affected at all by the electron-electron 
interaction within the wire |^ Q , and the shot noise of 
a wire with an impurity does not reflect the fractional 
charge of the backscattered quasiparticles Bil. The 
reason lies in the coupling of the QW to the FL leads 
at its ends: experimentally only the chemical potential 
of electrons in the leads can be manipulated, while the 
chemical potential of the fractionally charged quasiparti- 
cles is not controllable. This represents the main differ- 
ence with respect to the case of edge currents in fractional 
quantum Hall bars, where right and left movers are spa- 
tially separated (chiral LL), and the chemical potential 
of the Laughlin quasiparticles can be controlled by the 
upstream reservoir ^lOj . In contrast, in a semiconductor 
QW or a nanotube right and left movers interact and 
interfere in the same channel resulting in a nontrivial in- 
fluence of the contacts even for the case of ideal adiabatic 
coupling. 

Due to the important effects of the leads, a finite length 
QW in many respects acts as a complicated scatterer in 
a FL wire, and transport properties obey FL behavior. 
Thus, experimental observations of LL behavior in semi- 
conductor QWs [Til and nanotubes 13] a-re rather 
based on the temperature dependence of the conduc- 
tance, or on modifications by the Coulomb interaction 
of resonant tunnelling through quantum dot structures 



jlJI . Here we show that due to the interaction within the 
wire, the nonlinear conductance of a QW with a single 
impurity shows oscillatory behavior. In contrast to oscil- 
lations in double impurity systems ^Jl or in wires with 
tunneling contacts to the leads Il6j. the oscillations 
predicted here are completely absent in a non-interacting 
wire. So we do not discuss modifications of an oscillatory 
effect already present in FL theory, but non-monotonic 
behavior of the conductance that is a signature of LL 
behavior. 

To study a QW of finite length coupled to non- 
interacting reservoirs, we employ the inhomogeneous LL 
model I3 , where the interaction parameter has a spa- 
tial dependence. It varies from g{\x\ > L/2) = 1 in the 
leads to g{\x\ < L/2) = 5 < 1 in the QW. Although in 
real physical systems without backscattering at the con- 
tacts the variation ofg(x) from unity to g must be smooth 
on a length scale 3> Xp, where Af is the Fermi wave- 
length, it turns out that, as long as Lg ^ L, the simpler 
model of a g{x) step-function already captures the impor- 
tant physical properties of a QW adiabatically connected 
to non-interacting electron reservoirs. Moreover, we have 
to keep in mind that only under the assumption Lg <C i 
the QW has a well-defined length L. Throughout this 
article we treat the spin-polarized situation at zero tem- 
perature, leaving finite temperature aspects for a future 
publication. 

For fi=l, the full Lagrangian of our model is £ = £0 + 
Cj + Cij. In bosonization, see Ref. 5] for a recent review, 
the inhomogeneous LL Lagrangian of a clean system in 
terms of the phase field 6{x,t) reads 
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where wj? is the bare Fermi velocity of non- interacting 
right- and left-movers. The backscattering of an impurity 
at position xq yields the term 



Ci = —A cos 



4:Tr9{xo, t) + 2kFXQ 
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with the 'bare' impurity strength A, whereas the impu- 
rity forward scattering omitted here does not affect the 
current-voltage characteristics. A third term 



dxU{x)dxO , 



(3) 



describes the coupling to an electrostatic potential 
U{x), with eU{x < ~L/2) = ^l, eU{x > L/2) = ^ir, 
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and U{x) = for |a;| < L/2. Here, /Ul and /z_r are 
respectively the chemical potentials in the left and right 
reservoir, and we assume a symmetrically biased system 
with — eU/2 and = —eU/2. U{x) describes the 
shift of the band bottom due to electroneutrality in the 
FL leads 9J. In contrast, the QW itself does not remain 
electroneutral in presence of an applied voltage U and 
the electrostatics within an interacting QW emerges 
naturally as a steady-state effect 0, . 

In bosonization, the average current is given by (/) = 
— {e/y/TT)dt{9{x,t)), and it is most convenient to decom- 
pose the field operator 0{x,t) as |0 9{x,t) = 9p{x,t) + 
(j){x,t), where 9p{x,t) is a particular solution of the 
equation of motion following from the action S[9] = 
Sq + Si + Sij that corresponds to Eqs. 1^) to Q. On 
the other hand, t) can be regarded as the fluctua- 
tion operator which therefore has to fulfill 



dt{<l>{x,t))^0. 
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Straightforward algebra shows that the particular solu- 
tion may be written as 
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Consequently, the average current reads 

{I)-f{U-V), 

2,71 
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where is a parameter corresponding to the four- 
terminal voltage drop at the impurity site. In order to 
determine V , one has to impose the condition Q defining 
the fluctuation fleld ch. 



The effective action for is easily found to read 
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A / dt cos V4:T:4>{xo,t) - e{U -V)t + 2kFXQ 



For weak backscattering the terms related to the im- 
purity, namely the second and third line of Eq. l |7| l- 
can be treated as perturbations on the free boson part 
(first line) employing standard techniques. The con- 
ditions when such a perturbative approach is reliable 
will be discussed below. It is convenient to character- 
ize the conductance in terms of the backscattered cur- 
rent Ibs = Io- {!) = (eV27r)F, where Jo = {e^/2TT)U. 
As can easily be shown, the lowest order contribution 
to Ibs is of order A^. Without loss of generality we as- 
sume U > and find, using the requirement that the 
backscattered current takes the form 



-^BS — 



eX^gL 
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where u = eU/ujL is the ratio of the voltage to the fre- 
quency lul — VF/gL characterizing a ballistic traversal 
of the wire of length L. Here, the correlation function 
C{^:Vi''') — {4'{£,:T)4'{Vt^) ~ 0^('7jO))o in the absence of 
an impurity is expressed in terms of the dimensionless 
variables ^ = x/L, rj = y/L and r = tvF/gL- At the 
dimensionless impurity site one can show that 



where 7 = (1 — g)/{l + g), and a = u}l/uJc- Here luc is 
the standard high-energy cutoff of the LL model with A, 
ell, ujl Re- 



using Eqs. © and ®, the backscattered current 
can now be evaluated numerically for arbitrary values 
of the parameters g (interaction strength), u (voltage 
times length), and ^0 (impurity position). Fig. ^ depicts 
the result as a function of u for two impurity positions 
(in the middle and at 1/4 of L), and for an interaction 
value g — 0.25 characteristic for single-walled metallic 
carbon nanotubes [iS]. Apparently, the backscattered 
current oscillates as a function of u = eU/ujL, and the 
oscillations die out as u increases. 



Analytical expressions for the current can be obtained 
in special situations. For a very long wire (u — *■ 00) one 
can perform an asymptotic expansion yielding 



^BS =/|*s(C/)[l + fBi{u-^o) 



(10) 
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where /gg is the leading order term while /] 
the dominant oscillating correction. The dots represent 
further subleading orders. Explicitly one finds 



^BS(C^) 



T(2g) 27r (eC//wc)2(i-s) ' 
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which is independent of the length L of the wire and the 
impurity position, and exhibits a power-law behavior as 
a function of the applied voltage U in accordance with 
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the result for the homogeneous LL (see e.g. U). Af{a;g) 
is a constant of order unity. In contrast, the corrections 
exphcitly depend on both and L (through u), clearly 



J 



showing that oscillations are a finite-size effect. For < 
1^0 1 < 1/2 the asymptotic expansion yields 
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where !?((?; s^o) is a numerical factor of order unity. 
A similar asymptotic expansion has been performed in 
Ref. T^, where the crossover from LL to FL behavior 
was studied for a tunnel barrier. 

The two frequencies appearing in the cosine terms of 
Eq. (|12|) are related to the distances of the impurity from 
the two contacts (in units of the wire length L). The 
phase shifts Us — TTg{l ± 7/2) as well as the decay fac- 
tors of the oscillations directly depend on the electron- 
electron interaction strength. Importantly, when the in- 
teraction in the wire is switched off {g — > 1), i.e., when 
one deals with a homogeneous system, the oscillations 
vanish since T{g^) ^ oo in this limit. Clearly, 7bs also 
vanishes when there is no back-scatterer, i.e. A — > 0. 

One can understand the origin of these current os- 
cillations in terms of interference effects of plasmon 
modes which are reflected both by the impurity and, in 
an Andreev-type process 1^], by the wire- lead contacts. 
When varying either the length L of the wire, or the 
voltage (which appears as a phase shift in the effective 
action {T))), one basically changes the interference con- 
ditions, which results in the oscillatory behavior of the 
current. The period of the oscillations is only indepen- 
dent of g if the backscattered current is considered as 
a function of the dimcnsionless variable u; however, one 
has to keep in mind that u {— eUgL/vp) has an intrinsic 
vp/ g-dependence. 

We also point out that the expression H12(l is singular 
for ^0 ~* 0; this is just due to a mathematical compli- 
cation when the distances of the impurity from the two 
contacts become equal. Then pairs of poles in the cor- 
relator merge. The full expression ((HJ for the current is 
however perfectly regular, and indeed one can still calcu- 
late the asymptotic expansion for this particular impurity 
position, obtaining 
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It should be emphasized that in the asymptotic expan- 
sion (fTIHl the further subleading orders (represented by 
dots) can be shown to become increasingly important 
the stronger the interaction is, i.e., the lower the value of 
g. In practice, retaining only /gg and f^^g is a very good 
approximation for g > 0.5 and u ^ 10. Nevertheless, 
also for g < 0.5 and/or lower values of u (not too close 
to u = 1), Eq. l(Tn|l gives a qualitatively correct picture 
of the properties of the backscattered current as we have 
checked numerically. 
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FIG. 1: The backscattered current Ibs is depicted, in units 
of e{Xujl/uJ^)'^ /LOLJV{a; g), as a function of w = eU /ujl- The 
interaction strength is g = 0.25. The solid line corresponds to 
^0 = and the dashed line to ^0 = 0.25. Decaying oscillations 
with a constant period are clearly visible for the case ^0 = 0. 



We now discuss how the behavior of the backscattered 
current changes with the interaction strength g. In Fig. [21 
the full numerical result of Jbs , Eq. ^ , is plotted in units 
of /gg, Eq. (|11|1 . for three different values of g. For sim- 
plicity we have considered the case of an impurity located 
in the middle. The figure illustrates that the stronger 
the electron-electron interaction, the more pronounced 
the oscillations are. This is due to the fact that the 
Andreev-type reflections at the contacts are interaction- 
dependent with a reflection coefficient 7 = {\—g)/{\+g) 
[^. Thus, interference effects are most pronounced in 
the limit 5 — > 0. Furthermore, we see that the period 
of the oscillations is the same for any value of namely 
Au — 27r for ^0 = 0, but different values of g yield differ- 
ent phase shifts. 

In the opposite limit of a short wire (u — > 0), the in- 
tegral ||SJ) can also be evaluated analytically, and we find 
for the backscattered current 



lB5^n{X,g) — U 
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with the cut-off and interaction dependent 'reflection co- 
efficient' 



7^(A,g) =7r2AA2(a;g) 
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where the constant Af2{a; g) is of order unity. The power 
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FIG. 2: The backscattered current Ibs is plotted, in units 
of /gg, as a function of u = cU/lul for three different values 
of the interaction strength: g = 0.25 (dotted line), g = 0.50 
(dashed line), and g = 0.75 (solid line). Electron-electron 
interaction affects the decay as well as the phase shift of the 
oscillations. 



law behavior of 7?.(A, g) as a function of the system size L 
agrees with the calculation of the hnear conductance from 
the Kubo formula Furthermore, Ea. (|13|l bears a 

resemblance to the usual Landauer formula. In the limit 
g ^ 1 the reflection coefficient approaches the expected 
value 7^(A, 1) = n'^iX/iVc)^ . 

Finally we discuss the applicability of the perturbative 
approach. To be consistent, the backscattered current 
Ibs has to be small compared to the current Iq in the 
absence of an impurity, i.e., 

|/bs| «e2f//27r. (15) 

The analytical expressions obtained in the limits of both 



large u, Eqs. iCOl-OIJ' ^^'^ small u, Eqs. (QSJ-imi), sug- 
gest that the consistency condition can be expressed in 
a simple and cut-off independent way, if we introduce 
an effective impurity strength A* = LUc{X/iL!c)^''^^~^\ In- 
deed, Eq. ((T5|l then amounts to A* <^ eU for large u, 
and A* <C for small u, respectively. This is in ac- 
cordance with general reasoning. It is well known that 
the backscattering term Q is a relevant perturbation 
to the LL fixed point, and therefore the RG-flow would 
drive the system to two disconnected QWs, unless some 
energy-scale cuts off the flow of A. In our case the applied 
voltage U or the wire ballistic frequency lul can play such 
a role. Therefore, as long as A* is much smaller than ei- 
ther ell or ujl, the QW remains in the neighborhood of 
the LL fixed point and our perturbative treatment is jus- 
tified. In the limit 5 — > 1, the effective impurity strength 
A* vanishes, and perturbation theory is always applicable 
if A <C We- 
ill summary, we have investigated the backscattered 
current of an interacting QW of finite size L that is 
coupled to non-interacting electron reservoirs. We have 
found that the interplay between electron-electron inter- 
action and finite size gives rise to distinct features in the 
current-voltage characteristics. In particular, we have 
found an oscillation of the backscattered current as a 
function oi u = eU /ujl- For carbon nanotubes one can 
estimate a typical ballistic frequency of the order of a 
few meV. From Figs. Handle one can see that the oscil- 
lations are visible for eU ^ 10 ... 50 meV, a reasonable 
value, which is consistent with the cut-off of about 1 eV. 
Hence, the predicted effects should be observable. 
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